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FOREWORD 
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SUMMARY 


The  principal  objectives  of  the  current  project  included; 

1 .  Work  on  Torpedo  acetylcholinesterase  (TcAChE)  involving  two  principal  topics: 

A)  Preparation  and  characterization  of  'aged'  and  'non-aged'  organophosphoryl  (OP) 
conjugates. 

B)  Use  of  time-resolved  crystallography  in  order  to  understand  the  mechanism  of  action  of 
acetylcholinesterase  (AChE),  as  well  as  the  traffic  of  substrate  and  products  through  the 
active-site  gorge. 

2.  Structural  studies  on  other  cholinesterases  (ChEs). 

3.  Homology  model  building  and  detailed  stmctural  comparison  of  human,  bovine  fetal  semm, 
snake  venom  and  insect  AChE,  and  of  horse  and  human  serum  butyrylcholinesterase  (BChE). 

4.  Structural  studies  on  human  paraoxonase. 

In  the  first  year  of  this  project,  we  have  made  progress  towards  all  four  objectives,  as  well  as  in 

additional  areas  related  to  the  long-term  aims  of  our  research  program,  as  follows: 

A.  The  structures  of  aged  conjugates  obtained  by  reaction  of  TcAChE  with 
diwopropylphosphorofluoridate  (DFP),  O-wopropylmethylphosponofluoridate  (sarin),  and 
O-pinacolylmethylphosphonofluoridate  (soman)  were  solved  by  X-ray  crystallography  to 
2.3,  2.5  and  2.2k.  resolution,  respectively. 

B.  Laue  X-ray  data  were  collected  at  the  Third  Generation  Synchrotron,  i.e.  the  European 
Synchrotron  Research  Facility  (ESRF),  Grenoble,  on  trigonal  crystals  of  TcAChE  complexed 
with  edrophonium  (EDR).  Data  were  collected  in  msec  and  ps  time  slices;  even  with  such 
exposure  times,  it  was  possible  to  obtain  electron  density  maps  comparable  in  quality  to  those 
obtained  from  monochromatic  data. 

C.  The  3D  structure  of  a  complex  of  human  recombinant  AChE  with  the  snake  venom  toxin, 
fasciculin,  has  been  determined  to  2.1  k  resolution.  A  similar  data  set  has  been  collected  for  a 
complex  with  fasciculin  of  the  human  recombinant  E202Q  mutant. 

D.  Homology  models  have  been  constmcted,  on  the  basis  of  the  TcACtiE  structure,  for  human, 
Bungarus  fasciatus  and  Drosophila  melanogaster  AChE’s,  and  for  human  BChE.  Modelling 
reveals  that  the  overall  charge  asymmetry,  first  observed  for  the  Torpedo  enzyme,  is  conserved 
throughout  the  ChE  family. 

E.  A  number  of  neural  cell-adhesion  proteins,  including  gliotactin,  neurotactin,  and  neuroligin,  have 
substantial  sequence  homology  to  the  ChEs.  Homology  models  of  these  proteins  showed  that 
their  electrostatic  characteristics  were  similar  to  those  of  the  AChEs,  with  which  they  share  a 
negative  electrostatic  motif,  which,  in  AChE,  is  located  around  the  entrance  to  the  active-site 
gorge. 

F.  Cloning  and  expression  studies  of  human  paraoxonase  have  been  initiated. 
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OBJECTIVES 

The  principal  objectives  of  this  contract  are  to  better  understand  the  mechanism  of  action  of 
AChE,  and  to  characterize  its  various  ligand-binding  sites.  These  objectives  are  being  approached 
by  the  method  of  single  crystal  X-ray  diffraction,  together  with  time-resolved  X-ray  crystallography, 
using  Torpedo  califomica  AChE  (TcAChE)  as  the  principal  source  of  enzyme.  Computational 
molecular  model  building  and  electrostatic  characterization  are  being  used  as  complementary 
theoretical  approaches.  Attempts  are  also  being  made  to  crystallize,  and  if  successful,  to  solve  the 
3D  structures  of  cholinesterases  (ChEs)  from  other  sources.  These  include  human  recombinant 
AChE,  bovine  fetal  semm  AChE,  AChE  from  the  venom  of  the  kxait,  Bungarus  fasciatus  ffi/AChE), 
and  horse  serum  BChE.  Homology  model  building  will  be  used  to  compare  and  analyze  the 
stmctures  of  these  enzymes.  In  parallel,  we  are  attempting  to  express  and  crystallize,  with  the 
eventual  objective  of  solving  its  3D  stmcture,  human  semm  paraoxonase/arylesterase. 

1 .  Work  on  TcACtiE  involves  two  principal  topics: 

A)  Preparation  and  characterization  of  'aged'  and  'non-aged'  OP  conjugates,  with  the  objective 
of  fully  understanding  the  mechanistic  basis  and  specificity  of  interaction  of  organophosphates  and 
organophosphonates  of  therapeutic  and  toxicological  interest  with  the  active  site  of  AChE,  as  well  as 
the  molecular  basis  of  the  'aging'  process. 

B)  Use  of  time-resolved  crystallography  in  order  to  understand  the  mechanism  of  action  of 
AChE,  as  well  as  the  traffic  of  substrate  and  products  through  the  active-site  gorge.  The 
experimental  approaches  to  be  adopted  will  be  Laue  crystallography  and  ultra-rapid  freezing  time- 
resolved  crystallography.  In  either  case,  time-resolved  data  collection  will  be  preceded  by  synthesis 
and  characterization  of  suitable  caged  compounds,  and  by  a  study  of  their  interaction  with  TcAChE, 
initially  in  solution  and,  subsequently,  under  steady-state  X-ray  data  collection  conditions. 

2.  Work  on  other  ChEs  obviously  depends  on  the  success  of  crystallization  attempts.  Particular 
attention  is  being  paid  to  human  AChE  (liAChE),  in  view  of  its  toxicological  and  therapeutic 
importance.  In  the  first  year  of  the  contract,  we  have  devoted  substantial  efforts  to  improving 
purification  techniques  already  developed  in  our  own  laboratory  for  processing  culture  media 
containing  recombinant  hAChE  (rhACtiE),  supplied  by  the  group  of  Dr.  Avigdor  Shafferman  at  the 
Israel  Institute  for  Biological  Research  (IIBR). 

3.  Homology  model  building  and  detailed  stmctural  comparison  of  human,  bovine  fetal  semm, 
snake  venom  and  insect  AChE,  and  of  horse  and  human  semm  BChE. 

4.  Work  on  paraoxonase  is  subject  to  the  same  restrictions  as  for  the  various  ChEs,  namely  the 
availability  of  crystals.  Thus  expression  and  purification  will  precede  data  collection  and  stmcture 
determination. 
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BACKGROUND 


Acetylcholinesterase 

Acetylcholinesterase  (AChE)  plays  a  key  role  in  cholinergic  neurotransmission.  By  rapid 
hydrolysis  of  the  neurotransmitter,  acetylcholine  (ACh),  it  terminates  the  chemical  impulse,  thereby 
permitting  rapid  repetitive  responses  (1).  AChE  is,  accordingly,  characterized  by  a  remarkably  high 
specificity,  especially  for  a  serine  hydrolase  (2),  functioning  at  a  rate  approaching  that  of  a 
diffusion-controlled  reaction  (3, 4).  The  powerful  acute  toxicity  of  organophosphoms  (OP)  poisons 
(as  well  as  of  carbamates  and  sulfonyl  halides  which  function  analogously)  is  attributed  primarily  to 
the  fact  that  they  are  potent  inhibitors  of  AChE  (5,  6).  Inhibition  is  achieved  by  their  covalent 
attachment  to  a  serine  residue  within  the  active  site  (2,  7).  AChE  inhibitors  are  utilized  in  treatment 
of  various  disorders  such  as  myasthenia  gravis  and  glaucoma  (6,  8).  More  recently,  they  have  been 
under  active  consideration  for  use  in  the  management  of  Alzheimer's  disease  (9,  10).  The 
cholinesterase  inhibitors,  tacrine,  under  the  trade  name  Cognex®,  and  E2020,  under  the  trade  name 
Aricept®,  have  both  been  approved  for  use  by  the  FDA  (1 1,  12).  Elucidation  of  the  3D  structure  of 
AChE  is,  therefore,  of  fundamental  interest  for  understanding  its  remarkable  catalytic  efficacy,  and 
of  applied  importance  in  drug  design  and  in  developing  therapeutic  approaches  to  OP  poisoning. 
Furthermore,  information  concerning  the  ACh-binding  site  of  AChE  is  pertinent  to  an 
understanding  of  the  molecular  basis  for  the  recognition  of  ACh  by  other  ACh-binding  proteins 
such  as  the  various  ACh  receptors  (13-15). 

Crystallization  and  determination  of  the  3D  stmcture  of  AChE  from  TcAChE  (16,  17)  enabled 
us  to  visualize,  for  the  first  time,  at  atomic  resolution,  a  binding  pocket  for  ACh.  It  also  allowed  us 
to  identify  the  active  site  of  AChE,  which  is  located  at  the  bottom  of  a  deep  gorge  lined  largely  by 
aromatic  residues  (17,  18).  This  unusual  and  unexpected  structure  gives  us  the  opportunity  to  work 
out  structure-function  relationships  of  AChE.  Furthermore,  the  so-called  ‘anionic’  binding  site  for 
the  quaternary  moiety  of  ACh  does  not  contain  several  negative  charges,  as  was  earlier  postulated 
on  the  basis  of  the  ionic-strength  dependence  of  catalytic  activity  (19).  Our  modeling  studies  (17) 
suggested  that  the  quaternary  group  interacts  primarily  with  the  indole  ring  of  the  conserved 
tryptophan  residue,  W84;  this  was  borne  out  by  crystallographic  studies  on  several  AChE-ligand 
complexes,  complemented  by  photoaffinity  labeling,  which  also  suggested  a  prominent  role  for 
another  conserved  aromatic  residue,  F330  (20, 21). 

The  availability  of  the  3D  structure  of  7c AChE  stimulated  a  large  body  of  work,  in  several 
laboratories,  which  shared  the  common  objective  of  correlating  the  unexpected  structure  revealed  by 
X-ray  crystallography  with  the  catalytic  function  of  the  enzyme.  These  studies  can  be  loosely 
classed  in  four  main  categories,  two  theoretical  and  two  experimental  (some  studies  encompass 
more  than  one  category): 

1)  Theoretical  studies,  in  which  the  stmcture  of  AChE  or  BChE  from  various  sources  is 
modeled  on  the  basis  of  the  3-D  stmcture  of  TcAChE,  together  with  the  amino  acid  sequence  of 
the  ChE  in  question,  obtained  from  cloning  and/or  sequencing  studies.  Such  modeling  is 
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possible  due  to  the  high  sequence  homology  that  exists  between  ChEs  of  differing  phylogenetic 
origin  (22, 23). 

2)  Theoretical  studies  aimed  at  elucidating  the  electrostatic  and  dynamic  properties  of  the  ChE 
molecule  and  their  contribution  to  known  catalytic  properties. 

3)  Experimental  investigations  concerned  with  understanding  the  catalytic  and  ligand-binding 
properties  of  the  ChEs  on  the  basis  of  kinetic  and  spectroscopic  studies  with  AChE  and  BChE 
from  various  sources. 

4)  Studies  using  site-directed  mutagenesis  to  understand  the  contribution  of  particular  amino 
acid  residues  to  catalytic  activity,  substrate  specificity,  and  mode  of  interaction  with  reversible 
and  covalent  inhibitors. 

Inspection  of  the  3D  structures  of  TcAChE  and  of  various  AChE-ligand  complexes,  as  well  as 
of  the  stmcture  of  homologous  enzymes,  led  to  predictions  concerning  the  involvement  of  certain 
amino  acid  residues,  almost  all  localized  in  the  ‘aromatic  gorge’,  in  various  aspects  of  catalytic 
activity.  Site-directed  mutagenesis  was  the  principle  tool  for  testing  such  predictions.  Thus,  on  the 
basis  of  the  topology  of  the  gorge  (17),  of  the  X-ray  stmcture  of  complexes  of  TcAChE  with  the 
bisquatemary  inhibitors  decamethonium  (20)  and  BW284c51  (Harel  et  al,  unpublished),  and  on 
the  modeling  of  hBChE  (23),  it  was  suggested  that  the  so-called  ‘peripheral’  anionic  site, 
previously  identified  from  kinetic  (24)  and  spectroscopic  studies  (25,  26),  is  located  at  the  top  of  the 
gorge.  Like  the  ‘anionic’  subsite  of  the  active  site,  it  was  thought  to  contain  three  aromatic  residues, 
Y70,  Y121,  and  W279  {Torpedo  numbering  is  used  here  and  subsequently).  ICinetic  comparisons 
with  hBChE,  which  lacks  these  three  residues,  and  with  chicken  AChE,  which  lacks  two  of  them, 
provided  strong  experimental  evidence  to  support  this  prediction  (23,  27-30).  In  addition,  a  role  for 
the  acidic  residue,  D72,  in  the  ‘peripheral’  site,  was  proposed  on  the  basis  of  site-directed 
mutagenesis  (31).  Modeling  also  suggested  that  two  aromatic  residues,  F288  and  F290,  confer 
narrow  substrate  specificity  upon  AChE  compared  to  BChE  (32,  33);  this  suggestion,  too,  was 
clearly  borne  out  by  mutagenesis  studies.  Thus,  our  computer  modeling  had  high  predictive 
capacity,  and  we  are  now  modeling  hAChE  similarly,  as  well  as  more  distantly  related  enzymes,  e.g. 
insect  AChEs,  in  the  context  of  insecticide  development. 

Studies  carried  out  to  date  have  involved  about  200  mutations,  principally  of  hAChE,  mouse 
AChE  (mAChE),  TcAChE,  Torpedo  marmorata  AChE  (TmAChE)  and  Drosophila  melanogaster 
AChE  (DmAChE),  as  well  as  of  hBChE  (34).  An  updated  compilation  of  these  mutations  can  be 
obtained  by  accessing  the  Cholinesterase  Data  Base,  at  Montpellier,  France,  maintained  by  A. 
Chatonnet  (see:  http://www.ensam.inra.fr/cholinesterase/).  These  mutagenesis  studies  leave  open 
numerous  questions.  Thus,  the  present  data  ascribe  firm  roles  to  about  half  of  the  conserved 
aromatic  residues  in  the  gorge,  and  the  role  of  the  remainder  is  either  not  fully  established  or 
unknown  Furthermore,  the  role  of  two  acidic  residues,  E199  and  D72,  the  former  near  the  bottom 
of  the  gorge  and  the  latter  close  to  the  top,  is  not  clearly  ascertained,  and  the  long-standing  question 
of  the  molecular  basis  for  substrate  inhibition  is  still  unresolved.  The  attractive  possibility  that 
substrate  inhibition  involves  binding  of  a  second  ACh  molecule  at  the  ‘peripheral’  site  (27)  is 
prejudiced  by  a  number  of  studies.  These  involve  mutagenesis  of  aromatic  residues  at  the 
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'peripheral'  site  (28, 35)  and  comparative  studies  on  the  chicken  and  Torpedo  enzymes  (30),  as  well 
as  the  observation  that  mutation  of  El 99  strongly  reduces  substrate  inhibition  (36).  A  complication 
in  analyzing  all  the  site-directed  mutagenesis  data  is  that,  in  certain  cases,  mutations  at  the  top  of  the 
gorge  appear  to  affect  parameters  ascribed  to  structural  factors  associated  with  the  bottom  of  the 
gorge,  and  vice  versa  (31).  This  is  of  interest,  because  it  provides  evidence  that  the  gorge  may 
function  in  an  integrated  fashion,  but  also  means  that  interpretation  of  individual  site-directed 
mutagenesis  experiments  must  be  treated  with  caution. 

Paraoxonase 

Paraoxonase  (EC  3. 1.8.1)  belongs  to  the  family  of  A-esterases  as  defined  by  Aldridge  (37), 
which  are  capable  of  hydrolyzing  OPs,  carbamates  and  aromatic  carboxylic  acid  esters.  They  have 
been  the  subject  of  a  substantial  research  effort  due  to  their  potential  capacity  to  offer  protection 
against  intoxication  by  paraoxon  and  by  other  OPs,  which  might  serve  as  insecticides.  Indeed,  there 
is  evidence  to  support  the  notion  that  serum  paraoxonase  plays  a  significant  role  in  the 
detoxification  of  paraoxon  in  the  rabbit  in  vivo  (38).  Paraoxonase  may  also  confer  protection  on 
coronary  artery  disease  by  destroying  pro-inflammatory  oxidized  lipids  present  in  oxidized  low- 
density  lipoproteins  (LDLs).  Indeed,  paraoxonase  knock-out  mice  are  susceptible  both  to  OP 
toxicity  and  to  atherosclerosis  (39). 

OP-hydrolyzing  activity  has  been  described  throughout  the  phylogenetic  tree  from 
microorganisms  up  to  several  mammalian  species  (40,  41).  It  remains  to  be  established  whether 
various  esterases  capable  of  hydrolyzing  OP  in  different  organisms  e.g.  the 
diisopropylphosphofluoridate  hydrolase  (DFPase)  of  the  squid  head  ganglion  (42),  are 
homologous  members  of  a  small  number  of  families  displaying  substantial  structural  similarity.  In 
the  case  of  a  phosphotriesterase  from  Pseudomonas,  which  is  known  to  hydrolyze  both  paraoxon 
and  parathion  (43),  it  has  been  clearly  established  that  it  shows  no  homology  with  the  human 
enzyme,  which  does  not  hydrolyze  parathion  (44).  The  molecular  structure  of  the 
phosphotriesterase  has  recently  been  solved  to  2.1  A  resolution  (45,  46).  From  the  3D  structure, 
together  with  kinetic  and  chemical  data  on  this  enzyme  (47,  48),  it  does  not  appear  to  be  closely 
related  to  human  paraoxonase. 

The  human  serum  enzyme  which  hydrolyses  paraoxon  is  known  as  human 
paraoxonase/arylesterase  (PONl),  due  to  its  high  activity  on  aryl  esters  (49).  It  has  been 
investigated  extensively  by  the  laboratories  of  Furlong  (50)  and  La  Du  (44).  Thus  it  has  been 
purified  to  homogeneity  (51),  and  also  cloned  and  expressed  (52,  53).  It  occurs  in  two  isozymic 
forms,  A  and  B,  whose  pharmacogenetics  have  been  worked  out  by  the  two  laboratories.  The  two 
forms  differ  substantially  in  their  specifie  activity  towards  various  substrates.  Notably,  the  B-type 
hydrolyses  paraoxon  over  6-fold  faster  than  the  A-type  at  high  ionic  strength,  and  hydrolyses  2- 
naphthyl  acetate  2.5-fold  faster. 

The  laboratory  of  La  Du  has  developed  a  procedure  for  purification  of  PON  1  from  human 
serum  (51)  which  allowed  them  to  obtain  both  the  A-  and  B-type  enzymes  purified  over  500-fold 
(54).  PONl  binds  tightly  to  apolipoproteins,  and  also  displays  a  tendency  to  aggregate,  properties 
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which  necessitate  purification  in  the  presence  of  a  non-ionic  detergent  and  glycerol.  The  purified 
enzyme  is  stable  for  months  when  stored  at  4°C  in  the  presence  of  glycerol,  a  non-ionic  detergent, 
and  Ca”;  the  divalent  cation  is  essential  for  maintenance  of  the  native  conformation.  By  this 
procedure,  PONl  was  prepared  in  10  mg  batches,  which  permitted  preliminary  screening  of 
crystallization  conditions.  Indeed,  we  were  able  to  obtain  some  microcrystals  from  one  such  batch. 

The  polypeptide  chain  of  PONl  contains  354  amino  acids,  and  ca.  16%  carbohydrate,  which  is 
attached  at  three  glycosylation  sites.  La  Du  and  coworkers  (52)  detected  two  amino  acid  residues, 
54  and  191,  which  might  provide  the  structural  basis  for  the  A/B  polymorphism.  Subsequent 
molecular  genetic  and  nucleotide  sequencing  studies  by  the  La  Du  and  Furlong  groups  (53,  55) 
showed  that  the  Arg/Gln  polymorphism,  which  is  observed  at  residue  191,  follows  exactly  the  A/B 
phenotype  of  human  PON  1 .  No  influence  on  the  enzyme  has  yet  been  observed  of  the  Met/Leu 
polymorphism  which  is  commonly  observed  at  residue  54. 
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RESULTS 


Acetylcholinesterase  from  TcAChE 

a)  The  structure  of  'aged'  OP  conjugates  of  TcAChE 

OP  inhibitors  react  rapidly  with  AChE,  and  may  then  undergo  post-inhibitory  reactions  to 
produce  an  irreversibly  inhibited,  'aged'  OP-enzyme  conjugate.  To  understand  the  structural  basis 
for  the  non-reactivatibility  of  aged  OP-AChE  conjugates,  we  crystallized  and  solved  the  X-ray 
structures  of  conjugates  obtained  by  reaction  of  TcAChE  with  DFP,  sarin  or  soman.  The 
stoichiometric  conjugates  were  produced  in  collaboration  with  the  group  of  Dr.  Avigdor 
Shafferman  at  IBBR.  After  they  had  been  shown  to  be  fully  inactivated  (>99%)  and  aged  (>95%), 
each  conjugate  was  crystallized  from  PEG200  under  conditions  similar  to  those  employed  for 
native  TcAChE  (56). 

Two  different  crystal  forms  suitable  for  X-ray  diffraction  studies  appeared  in  drops  containing 
the  aged  OP- TcAChE  conjugates:  orthorhombic  space  group  P2,2,2|  and  trigonal  space  group 
P3,21.  X-ray  data  were  collected  using  trigonal  crystals,  and  finally  refined  at  2.4A  (DFP),  2.5 A 
(sarin),  and  2.2A  (soman)  resolution.  The  highest  positive  difference  density  peak,  corresponding  to 
the  OP,  was  observed  to  be  within  covalent  bonding  distance  of  the  Oy  of  Ser200  in  each  structure 
(Fig.  1).  The  final  reaction  product  was  methylphosphonylated  (MeP)-AChE  for  both  soman  and 
sarin  (Fig.  2),  and  monowcpropylphosphorylated  (MiPrP)-AChE  with  DFP  (Fig.  3). 

The  large  molar  excess  of  sarin  or  soman  used  should  have  resulted  in  almost  exclusive  reaction 
of  AChE  with  the  toxic,  P(S)  stereoisomers  (57).  Assuming  no  reorientation  of  the  inhibitor  in  the 
active  site  during  aging,  the  absolute  configuration  of  the  OP  adduct  in  the  MeP-AChE  structures  is 
consistent  with  a  simple  in-line  nucleophilic  attack  and  inversion  of  stereochemistry  at  the 
phosphorus  atom  during  phosphonylation.  The  inversion  mechanism  has  been  proposed  for  AChE 
by  analogy  with  other  serine  hydrolases  (58).  Comparison  of  the  aged  sarin  and  DFP  conjugates, 
both  of  which  have  lost  an  Aopropyl  group  during  aging,  supports  this  conclusion. 

The  primary  problem  created  by  the  dealkylation  reaction  of  aging  is  the  introduction  of  a 
formal  negative  charge,  because  anionic  phosphoesters  are  inherently  resistant  to  nucleophilic  attack 
(59).  Specific  interactions  between  the  enzyme's  active  site  and  the  OP  almost  certainly  add  to  the 
stability  of  aged  AChE,  however,  beeause  denaturation  of  aged  M/PrP-AChE  does  permit  base- 
catalyzed  dephosphorylation  (60). 

A  hydrogen  bond  from  Ne2  of  the  active  site  imidazolium  to  one  oxygen  atom  of  the  anionic 
OP  conjugate  has  been  postulated  to  be  a  key  interaction  that  stabilizes  both  aged  MiPrP-trypsin 
(61,  62)  and  M/PrP-chymotrypsin  (63).  Active  site  distances  and  angles  of  the  three  solved  OP- 
AChE  structures  are  consistent  with  the  formation  of  such  an  H-bond. 

Our  results  reinforce  the  conclusion  that  there  are  multiple  factors  militating  against  reactivation 
of  aged  AChE:  (i)  the  anionic  character  of  the  phosphoester;  (ii)  three  potential  H-bond  donors 
from  the  oxyanion  hole  to  one  of  the  oxygen  atoms  bound  to  P;  (iii)  an  additional  H-bond  from  the 
His440  N82  to  the  other  oxygen  atom;  and  (iv)  the  wrong  protonation  state  and  poor  geometry  for 
His440  to  act  as  a  general  base  catalyst. 
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The  phosphonyl  oxygen  atoms  may  cany  significantly  greater  fractional  negative  charges  than 
do  the  oxygens  of  the  ACh  transition  state  (TS)  adduct  {cf.,  (58)  and  (62)).  Despite  this  limitation, 
the  active  sites  of  MeP-AChE  should  provide  useful  structural  analogs  for  modeling  the  deacylation 
TS  formed  by  the  natural  substrate  (ACh)  (62, 64,  65). 

The  position  of  the  phosphonyl  oxygen  should  mimic  that  of  the  TS  carbonyl  oxygen.  All  three 
OP  structures  corroborate  that  the  oxyanion  is  stabilized  in  AChE  by  potential  H-bonds  from  the 
amide  nitrogens  of  G118,  G119  and  A201  (Figs.  2,  3).  This  is  consistent  with  the  preliminary 
model  for  the  binding  of  ACh  to  native  AChE  (17),  as  well  as  with  the  structure  of  a  complex  of  a 
transition  state  analog,  namely  m-(N,N,N-trimethylammonio)trifluoroacetophenone  (TMTFA)  with 
TcAChE  (21). 

A  hydrophobic  pocket  which  surrounds  the  phosphonate  methyl  group  is  formed  in  TcAChE 
by  W233,  F288,  and  F290.  F331  may  contribute  to  the  acyl  pocket  by  providing  7t-TC-interactions 
with  F288.  Both  F288  and  F290  provide  close  non-bonded  contacts  to  the  sarin  or  soman  methyl 
group  (Fig.  2),  suggesting  that  these  aromatic  side  chains  form  a  specific  pocket  for  the  methyl  of 
the  acetyl  group  of  ACh.  Contacts  from  the  acyl  pocket  to  the  phosphonate  methyl  group  are 
generally  closer  than  those  observed  for  the  CF3  group  of  TMTFA. 

Structure-activity  studies  with  carboxyl  ester  substrates  and  OP  inhibitors  first  suggested  the 
presence  of  specific  binding  pockets  in  the  active  site  of  AChE  that  complement  the  acyl  and  alkoxy 
groups  (reviewed  in  (66,  67)).  Furthermore,  this  approach  showed  convincingly  that  the  acyl  pocket 
of  AChE  was  more  restricted  than  that  of  the  closely  related  enzyme,  BChE.  When  the  crystal 
structure  of  AChE  became  available,  homology  models  of  BChE  were  built  and,  by  inference  from 
site-specific  mutagenesis  results,  two  aromatic  residues,  F288  and  F290,  were  implicated  as  the 
source  of  greater  selectivity  of  AChE  (29,  33,  36).  The  pocket  of  BChE  was  viewed  as  more  "open" 
because  it  contained  Leu  and  Val  (or  lie)  residues  at  these  positions,  whereas  the  two  Phe  residues 
in  the  pocket  of  AChE  were  described  as  a  "clamp"  for  the  methyl  group  of  ACh  that  would  poorly 
accommodate  bulkier  acyl  groups  (36). 

DFP  binds  >150-fold  more  tightly  to  BChE  than  to  AChE  (68).  Using  mutagenesis  and  a  series 
of  OPs,  it  was  shown  that  the  relatively  poor  binding  of  DFP  could  be  accounted  for  almost 
completely  by  an  increase  in  the  of  the  hAChE-DFP  reversible  complex.  This  increase  was 
mediated  primarily  by  F295  and,  to  a  lesser  extent  by  F297,  corresponding  to  residues  F288  and 
F290,  respectively,  in  TcAChE  (69). 

We  observed  the  first  major  conformational  change  in  a  TcAChE  complex  for  the  M/PrP-AChE 
structure  (obtained  after  reaction  with  DFP;  Fig.  3).  Significant  movement  of  the  main  chain  atoms 
allows  repositioning  of  F288  and  F290  to  accommodate  the  remaining  Aopropyl  group  of  the  OP 
in  the  'aged'  conjugate.  Superposing  the  MiPrP-AChE  structure  with  native  TcAChE  showed  that 
DFP  causes  the  Ca  atoms  at  positions  288-290  to  move  more  than  3.0A.  The  free  energy  penalty 
associated  with  this  movement  probably  explains  the  reduced  binding  affinity  of  AChE,  relative  to 
BChE,  for  DFP,  the  P(R)  phosphonates,  and  other  substrates  or  inhibitors  with  bulky  acyl  groups. 
The  conformational  change  in  the  acyl  pocket  observed  by  X-ray  crystallography  belies  earlier 
molecular  modeling  studies,  which  generally  predicted  that  the  degrees  of  freedom  available  to  the 
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substrate  or  inhibitor  were  restricted  by  F288  and  F290,  and  not  that  AChE  itself  undergoes  a  major 
conformational  movement  of  the  main  chain  to  accommodate  bulky  ligands. 

Although  there  is  evidence  for  conformational  changes  associated  with  aging  of  some  OP- 
AChE,  -BChE,  or  -chymotrypsin  conjugates  in  solution  (70-72),  we  found  no  major  structural 
rearrangements  in  TcAChE  as  a  result  of  the  aging  reaction  itself.  The  root  mean  square  deviation 
found  by  comparing  the  Ca  atoms  of  either  MeP-AChE  structure  with  native  AChE  was  less  than 
0.3A.  Moreover,  superposing  the  active  site  residues  of  each  MeP-AChE  with  those  of  native 
2ACE,  as  well  as  with  thirteen  other  previously  solved  TcAChE  stmctures,  showed  no  significant 
movements  that  could  be  attributed  to  aging.  It  remains  possible,  however,  that  a  transient 
conformational  change  occurs  on  the  pathway  to  the  final  aged  structure  observed  by  X-ray 
crystallography.  This  possibility  may  be  investigated  by  solution  studies,  as  well  as  by  ongoing  X- 
ray  crystallographic  studies  of  non-aged  OP- AChE  complexes. 
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b)  X-ray  Laue  studies  on  TcAChE 

As  already  mentioned,  AChE  is  an  extremely  efficient  catalyst,  with  a  turn-over  number  of 
"20,000  per  second  (2).  The  X-ray  structure  of  TcAChE  (17)  unexpectedly  revealed  that  the  active- 
site  is  buried  within  the  protein,  apparently  accessible  only  via  a  long  and  narrow  gorge  lined  with 
aromatic  residues.  Furthermore,  the  AChE  molecule  displays  an  asymmetric  charge  distribution, 
resulting  in  a  strong  dipole  moment,  aligned  approximately  along  the  active-site  gorge  (73,  74).  It 
was  suggested  that  the  electric  field  so  generated  attracts  the  positively  charged  ACh  to  the  rim  of 
the  gorge  (73).  The  substrate  is  then  guided  by  the  electric  field,  using  the  aromatic  residues  in  the 
gorge  as  a  series  of  low-affinity  sites,  towards  the  active  site.  The  narrow  dimensions  of  the  active- 
site  gorge,  taken  together  with  the  potential  gradient  along  it,  raise  cogent  questions  regarding  traffic 
of  substrate  and  products  to  and  from  the  active  site.  Thus,  in  addition  to  steric  considerations,  the 
potential  along  the  gorge  would  be  expected  to  accelerate  expulsion  of  acetate  from  the  active  site, 
but  retard  release  of  choline.  Furthermore,  movement  of  water  between  the  bulk  phase  and  the  gorge 
must  be  taken  into  account.  It  was,  therefore,  suggested  that  the  reaction  products,  choline  in 
particular,  might  leave  the  active  site  through  a  "back  door",  a  transient  opening  in  the  thin  aspect  of 
the  gorge  wall,  produced  by  movement  of  one  or  more  amino  acid  residues  (73,  75,  76).  Residues 
that  might  contribute  to  the  opening  of  such  a  putative  hack  door  would  most  likely  be  located  on 
the  Q  loop  between  C67  and  C94  in  TcAChE  (18),  a  loop  which  is  known  to  adopt  an  open  and 
closed  conformation  in  the  stmcturally  similar  enzyme,  Candida  rugosa  lipase  (77).  Further 
evidence  for  an  alternative  entrance  to  the  gorge  comes  from  studies  on  the  mode  of  inhibition  of 
AChE  by  the  potent  snake  venom  polypeptide  toxin,  fasciculin.  Crystallographic  studies 
demonstrate  that  the  toxin  completely  blocks  the  entrance  to  the  active-site  gorge  (78,  79). 
Nevertheless,  significant  residual  catalytic  activity  can  be  observed  in  solution  (80).  However,  site- 
directed  mutagenesis  studies  which  were  designed  to  prove  the  importance  of  the  mobility  of  certain 
residues,  did  not  provide  support  for  the  existence  of  a  back  door  (76,  81). 

Time-resolved  crystallography  is  clearly  an  approach  with  the  potential  to  provide  direct 
experimental  information  to  either  support  or  refute  the  existence  of  a  “back  door”.  The  Laue 
diffraction  method  has  been  employed  successfully  in  recent  years  to  observe  enzyme  reactions  as 
they  occur  in  the  crystal  state  (82,  83).  The  technique  is  attractive  due  to  the  very  short  exposure 
times  required.  Recently,  crystallographic  snapshots  of  structural  changes  occurring  in  a  protein  on 
a  time  scale  as  short  as  nanoseconds  have  been  reported  (84).  A  prerequisite  for  successful  Laue 
measurements  is  achievement  of  a  high  degree  of  synchronization  of  the  initialisation  of  the 
reaction.  In  some  cases,  the  substrate  is  immobilized  by  attachment  of  a  photosensitive  chemical 
moiety,  effectively  turning  it  into  an  inhibitor.  Such  caged  compounds  can  then  be  activated  by 
photolysis,  usually  with  laser  flashes.  In  the  case  of  AChE,  a  series  of  caged  compounds  has  been 
synthesized,  and  characterized  photochemically  and  enzymatically,  which  generate  either  choline  or 
carbamylcholine  upon  photolysis  (85,  86).  CCh  is  a  substrate  of  AChE;  the  first  step  of  the 
reaction,  carbamylation  of  the  enzyme,  is  rapid,  but  decarbamylation  is  slow,  on  a  time  scale  of 
minutes  at  room  temperature  (87).  It  is  envisaged  that  synchronized  photolysis  of  caged  CCh  within 
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the  crystal  will  produce  rapid  release  of  choline  by  the  enzymatic  machinery,  and  that  movement  of 
choline,  and  concomitant  changes  in  the  protein  conformation,  may  be  monitored  by  the  Laue 
method. 

The  Laue  technique  has  some  inherent  difficulties,  such  as  a  paucity  of  low-resolution 
reflections,  overlap  of  multiple  spots,  spatial  overlap  of  streaked  spots,  and  an  unfavorable  peak-to- 
background  ratio  (88).  These  drawbacks  result  in  a  requirement  for  well-diffracting  crystals  with 
low  mosaic  spread.  Since  success  in  Laue  diffraction  experiments  is  not  guaranteed  even  if  good 
diffraction  is  obtained  under  conditions  of  a  monochromatic  experiment,  evaluation  of  the  quality  of 
the  data  obtainable  by  the  Laue  technique  for  trigonal  crystals  of  TcAChE  was  undertaken.  The 
initial  goal  was  to  obtain  a  data  set  of  sufficient  quality  to  permit  resolution  of  structural  features 
within  the  active  site.  To  test  this,  crystals  were  soaked  with  a  solution  of  a  small  reversible  ligand, 
edrophonium  (EDR)  (see  Fig.  4).  The  stmcture  of  the  complex  of  TcAChE  with  EDR  was 
determined  earlier  using  monochromatic  radiation  (20).  Comparison  of  electron  density  maps 
obtained  from  Laue  data  with  maps  obtained  from  monochromatic  data  should  reveal  whether  use 
of  the  Laue  technique  is  feasible  for  performing  time-resolved  studies  on  AChE. 

The  monochromatic  data  set  on  the  TcAChE-EDR  complex  that  was  used  as  a  reference  for 
estimating  the  quality  of  the  Laue  data  set  was  obtained  from  crystals  grown  in  the  same  way  as 
those  for  the  Laue  data  set,  using  PEG200  as  a  precipitant  (56),  rather  than  the  ammonium  sulfate 
employed  previously  (20).  The  new  data  were  collected  at  4°C,  at  the  EMBL  outstation  of  the  DESY 
in  Hamburg,  Germany,  beam  line  XI 1,  on  a  30cm  MAR  Research  IP  detector  (X.=0.92A).  The 
crystal  had  a  mosaic  spread  of  0.29°  and  diffracted  to  a  resolution  of  2.4A;  further  details 
conceeming  the  quality  of  the  data  are  given  in  Table  1.  The  structure  obtained  from  the 
monochromatic  synchrotron  data  is  virtually  identical  to  the  previously  determined  structure  with 
PDB  access  code  lACK  (20);  the  r.m.s.d.  between  coordinates  for  all  Ca  atoms  is  0.29 A,  and  a 
plot  of  the  r.m.s.d.  values  of  the  backbone  versus  residue  number  does  not  show  any  major 
conformational  changes.  The  r.m.s.d.  for  the  ten  non-hydrogen  atoms  of  the  ligands  is  0.30A. 

Laue  data  on  crystals  of  TcAChE  soaked  with  EDR  were  collected  at  the  European  Synchrotron 
Radiation  Facility  (ESRF)  in  Grenoble,  France,  at  the  white  radiation  beam  line  BL3  (ID9),  with  a 
wavelength  range  of  0.4  to  1.5  A.  The  crystals,  mounted  in  a  capillary,  were  tested  at  4°C,  using  an 
FTS  cooling  device.  Several  crystals  had  to  be  discarded  due  to  lack  of  high-resolution  reflections 
and/or  severe  streaking  of  the  spots.  Good  diffraction  was  obtained  from  a  rather  large  crystal, 
0.5x0.5x0.4  mm^,  that  showed  a  reasonably  low  mosaicity.  24  still  frames  with  an  exposure  time  of 
1  ms  and  a  crystal-to-film  distance  of  300  mm  were  collected  while  rotating  the  crystal  over  5° 
along  the  spindle  axis  between  exposures.  The  data,  collected  on  an  X-ray  image  intensifier  coupled 
to  a  CCD  camera,  were  corrected,  using  the  package  FIT2D  (89),  for  spatial  distortions  and  non- 
uniformity  of  response. 

The  crystal  orientation  parameters  and  relative  cell  dimensions  were  determined  using  the 
program  LAUECELL  (90).  The  scaled  cell  dimensions  (a=b=113.0,  c=136.8A)  are  close  to  the 
values  found  for  monochromatic  data  on  non-frozen  trigonal  TcAChE  crystals  in  general.  The 
program  LAUE  VIEW  (91)  was  used  for  further  processing  of  the  data.  The  Laue  spots  on  the 
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patterns  were  predicted  using  a  resolution-dependent  wavelength  bandpass,  based  on  a  known 
starting  A,-curve  and  Wilson  statistics.  Since  the  distribution  of  the  single  spots  as  a  function  of  the 
wavelength  is  not  uniform,  wavelength  normalization  was  done  in  two  parts  (see  Fig.  5).  The  single 
spots  are  more  densely  distributed  between  0.45A  and  0.9A  than  at  wavelengths 

higher  than  0.9 A.  The  A.-curve  for  wavelengths  between  0.45  and  0.9 A  was  simulated  with  a  16- 

o 

term  Chebychev  polynomial,  and  the  curve  for  wavelengths  between  0.9  and  1.25 A  with  a  3-term 
Chebychev.  Fig.  5  also  shows  a  Bragg-angle  cut-off  (i.e.  no  spots  in  the  upper  right  corner  of  the 
figure):  the  crystal  diffracted  beyond  the  edge  of  the  detector.  However,  the  risk  of  overlapping 
spots  prevented  a  significant  reduction  of  the  crystal-to-detector  distance. 

o  o 

The  resolution  limit  (d^j„)  used  for  processing  the  Laue  data  set  is  2.8A,  compared  to  2.4A  for 
the  monochromatic  data.  The  combined  X,-curve  (shown  in  Fig.  6)  was  used  for  deconvolution  of 
the  multiples  (92);  this  significantly  increased  the  completeness  of  the  data,  especially  in  the  oo- 
2dn.i„  range. 

Details  of  the  quality  of  the  Laue  data  set  are  given  in  Table  2.  The  amplitudes  of  the  full  Laue 
data  set  give  a  correlation  of  91.6%  with  those  of  corresponding  reflections  in  the  monochromatic 
data  set  (see  Fig.  7);  for  the  singles,  the  correlation  is  93.4%  and  for  the  deconvoluted  multiples, 
87.8%. 

The  same  refinement  procedure  was  used  for  both  the  monochromatic  and  Laue  data  sets.  The 
coordinates  of  native  TcAChE  in  the  trigonal  crystal  form  (PDB  access  code  2ACE)  were  used  as  a 
starting  model;  5%  of  the  data  (1891  reflections  for  the  monochromatic  and  1032  for  the  Laue  data) 
were  reserved  for  use  in  a  test  set  for  cross-validation  of  the  refinement  process.  To  reduce  the 
chance  of  introducing  model  bias  into  the  initial  maps,  the  reflections  for  the  test  set  were  taken  by 
selecting  the  corresponding  reflections  from  the  test  set  used  for  refinement  of  the  native  structure 
(56).  Initial  maps  were  calculated  after  rigid-body  movement  of  the  coordinates  for  protein  residues 
using  the  program  X-PLOR  v.  3.851  (93).  The  difference  between  R  and  Rf^^  at  this  point  was 
2.6%  for  the  monochromatic  and  2.2%  for  the  Laue  data,  compared  to  3.6%  for  the  final  stage  of 
refinement  of  2ACE  without  waters.  A  random  test  set  would  have  given  R  and  of  the  same 
value;  the  obtained  differences  indicate  that  most  of  the  model  bias  has  been  eliminated.  Subsequent 
refinement  was  done  using  the  maximum-likelihood  refinement  procedure  of  REFMAC  (94),  using 
all  data  between  20.0A  and  the  highest  resolution,  and  manual  model  building  using  O  (95).  The  R- 
factor  (Rfree)  ^fter  refinement  is  21.3%  (25.8%)  for  the  monochromatic  and  20.7%  (26.8%)  for  the 
Laue  data.  The  final  model  contains  252  waters  for  the  monochromatic  structure  and  144  for  the 
Laue  stmcture,  and  the  r.m.s.  deviation  between  the  two  structures  is  0.20A  for  all  Ca  atoms. 

Obtaining  a  data  set  of  reasonable  quality  from  Laue  diffraction  is  not  a  guaranteed  success  for 
every  crystallized  protein.  The  Laue  method  has  been  applied  with  good  results  to  well-behaving 
proteins  such  as  cutinase  (96),  restrictocin  (97)  and  for  photoactive  yellow  protein  (98);  it  has  also 
been  applied  to  crystals  of  high  symmetry,  such  as  tomato  bushy  stunt  vims  (99).  Since  the 
question  of  the  exit  route  for  reaction  products  from  the  active  site  of  AChE  had  been  raised,  we  had 
sought  a  method  that  would  permit  us  to  establish  unequivocally  whether  or  not  there  is  a  "back 
door"  to  the  active-site  gorge.  Laue  diffraction,  with  its  extremely  short  exposure  times  and  high 
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information  content  per  image,  could  well  supply  such  information,  provided  that  four  conditions 
are  met: 

•  initiation  of  the  enzymatic  reaction  must  be  synchronised 

•  the  reaction  needs  to  be  on  a  time  scale  commensurate  with  that  of  data  collection 

•  the  existence  of  intermediate  states  needs  to  be  ascertained 

•  the  Laue  data  need  to  be  of  adequate  quality 

Synchronization  is  currently  being  investigated  using  caged  compounds  especially  developed 
and  synthesized  for  this  purpose.  The  compounds  synthesized  are  nitrobenzyl  derivatives  of  choline 
and  of  carbamylcholine  (CCh),  which  generate  the  corresponding  free  ligands  upon  photoactivation 
(85,  86).  Their  interaction  with  TcACtiE  has  been  studied  both  in  solution  (86),  and  in  the  crystal 
state,  using  monochromatic  crystallography  (Raves,  unpublished  results). 

The  existence  of  one  intermediate  state  in  the  hydrolysis  of  ACh,  the  acetylated  enzyme,  is  well 
known  (100).  The  carbamylated  enzyme  formed  during  the  hydrolysis  of  CCh  provides  a  more 
stable  intermediate  that  should  be  observable  using  time-resolved  crystallography,  because  the 
decarbamylation  of  the  enzyme  through  hydrolysis  is  much  slower  than  deacylation  (87).  However, 
since  the  pathway  of  release  of  the  reaction  product  of  interest,  choline,  is  unknown,  the  question  of 
whether  an  intermediate  state  is  involved  in  the  release  cannot  be  answered  at  this  point. 

The  quality  of  the  Laue  data  depends  on  the  resolution  and  completeness  attainable,  the 
rehability  of  the  measured  intensities,  and  the  interpretability  of  the  electron  density  maps  that 
follows  from  these  factors. 

Unlike  in  monochromatic  diffraction  experiments,  the  high-resolution  spots  in  a  Laue  image  not 
only  have  a  higher  Bragg  angle  on  average,  placing  them  further  towards  the  edge  of  the  detector, 
they  also  cause  a  higher  density  of  spots  in  all  regions  of  the  image.  In  the  case  of  AChE,  where  all 
three  unit  cell  axes  are  longer  than  lOOA,  and  the  mosaicity  of  the  crystals  is  rarely  below  0.2°, 
spatial  overlap  of  spots  is  problematic.  Although  the  program  used  for  integration,  LAUEVIEW 
(91),  is  able  to  deconvolve  overlapping  reflections  to  some  extent,  simulations  showed  that  a 
serious  reduction  of  the  crystal  to  detector  distance  (e.g.  200  instead  of  300  mm)  was  unlikely  to 
give  useful  data  for  these  crystals.  A  theta  offset  of  the  detector  could,  however,  have  been 
instmmental  in  obtaining  more  unique  data,  as  is  indicated  by  the  Bragg  cut-off  in  Fig.  5.  Another 
data  collection  strategy  that  could  have  been  beneficial  for  these  crystals  was  described  by  Yang  et 
al.  (97):  to  reduce  the  extent  of  X-ray  heating  of  the  crystal  by  regulating  the  exposure  time.  The 
resolution  cut-off  that  was  chosen  for  the  data  set  described  in  the  previous  section  was  2.8A.  For 
the  purpose  of  these  studies  and  the  kind  of  structural  changes  we  hope  to  see,  2.8A  resolution 
should  be  sufficient. 

In  terms  of  completeness  (see  Table  2)  the  quality  of  the  dataset  is  veiy  good:  95.2%  in  the 
medium  resolution  range  (5.6-3.4A).  The  problem  of  the  low-resolution  hole  (88)  has  largely  been 
overcome,  given  the  completeness  of  82.8%  between  oo  and  5.6A  (2d^|„),  which  stems  from  the  fact 
that  a  relatively  large  number  (24)  of  frames  were  collected,  aided  by  the  powerful  deconvolution 
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algorithm  of  the  program  that  was  used  for  processing  the  data,  LAUEYDEW  (91).  At  higher 
resolutions,  the  completeness  of  the  reflections  is  lower,  but  adequate  for  inclusion  in  the  dataset. 

The  amplitudes  of  the  reflections  correlate  fairly  well  with  those  of  the  monochromatic  data  set. 
From  Fig.  7  it  can  be  seen  that  there  is  a  slight  tendency  to  overestimate  weak  reflections,  but  the 
overall  correlation  coefficient  of  91.6%  indicates  that  the  amplitudes  are  fairly  reliable.  The  of 
the  Laue  data  set,  14.8%  vs.  6.5%  for  the  monochromatic  data  up  to  2.8A,  also  indicates  a  reduction 
of  the  quality  of  the  data;  but  is  still  acceptable  given  the  reduction  in  total  X-ray  exposure  from 
several  hours  to  24  ms.  Since  the  program  LAUEVIEW  (91)  gives  the  user  full  manual  control  of 
the  integration,  rejection  and  scaling  criteria,  the  completeness  of  the  dataset  is  a  trade-off  against 
other  parameters  which  are  indicative  for  the  quality  of  the  data.  The  inclusion  of  the  non-redundant 
single-observation  reflections  into  the  dataset,  for  example,  improved  the  completeness,  but  caused 
some  deterioration  in  the  correlation  with  the  monochromatic  data.  In  addition,  the  determination  of 
the  wavelength  normalisation  curve  strongly  depends  on  the  accuracy  of  the  intensities,  and  special 
care  had  to  be  taken  not  to  include  too  many  unreliable  measurements,  while  still  being  left  with  as 
high  a  redundancy  as  possible.  The  final  wavelength  normalisation  curve  obtained  (Fig.  6)  did  not 
show  details  such  as  the  Pt-edges  arising  from  the  focussing  mirror  (91),  but  has  the  correct  overall 
appearance. 

Special  consideration  was  given  to  the  problem  of  model  bias  in  the  initial  electron  density 
maps,  since  the  structure  to  be  determined  from  the  Laue  data  is  very  similar  to  that  of  the  protein  in 
the  native  state.  For  this  purpose,  the  same  reflections  that  were  used  in  the  test  set  in  the  refinement 
of  the  initial  starting  model  (56)  were  flagged  in  the  monochromatic  and  Laue  data  sets  of  the  EDR 
complex  (101).  After  rigid-body  refinement  of  the  entire  protein,  the  difference  between  R  and  R^ 
indicates  that  most  of  the  model  bias  had  been  avoided. 

Fig.  8  shows  a  comparison  of  the  Fourier  difference  maps,  after  rigid-body  refinement  of  the 
active-site  region,  for  the  Laue  data  set  vs.  the  complete  monochromatic  data  set,  and  vj'.  the 
selection  of  reflections  from  the  monochromatic  data  set  which  are  present  in  the  Laue  data.  As  is  to 
be  expected,  the  electron  density  difference  maps  made  using  the  Laue  data  are  of  somewhat  lower 
quality  than  the  ones  obtained  using  the  monochromatic  data;  the  signals  are  weaker,  and  the  maps 
are  noisier  and  show  less  detail,  reflecting  the  difference  in  R^^.^  values  and  completeness  for  the 
two  data  sets.  Nevertheless,  the  difference  map  obtained  from  the  Laue  data  clearly  shows  the 
location,  the  orientation  and  the  shape  of  the  soaked-in  inhibitor,  EDR. 

In  addition,  a  peak  search  on  the  initial  Laue  difference  map  yielded  78  peaks  higher  than  3.0  0 
that  correspond  to  water  molecules  in  the  monochromatic  stmcmre,  indicating  that  the  maps  are 
interpretable.  For  some  of  the  incomplete  side  chains  in  the  initial  model  (61  atoms  in  23  different 
residues  are  disordered  in  the  native  stmcture),  several  could  be  located  in  the  maps  for  the 
monochromatic  stmcture,  reducing  the  number  of  missing  atoms  to  57  in  22  residues.  The  Laue 
maps  showed  even  more  of  the  missing  atoms,  leaving  only  51  atoms  in  19  residues  undetermined. 
It  is  especially  encouraging  to  observe  that  the  conformational  change  that  Oy  of  the  active-site 
S200  undergoes  upon  binding  of  EDR  (20),  from  hydrogen-bonding  to  H440  in  the  native  state  to 
interacting  with  the  oxyanion  hole  in  the  complex,  is  visible  in  the  electron-density  maps  of  the  Laue 
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data  set  (see  Fig.  9).  Such  detailed  information  will  be  extremely  relevant  in  the  determination  of  a 
transient  intermediate  structure  and  of  structural  changes  during  the  enzymatic  reaction. 

Based  on  the  data  presented  in  this  chapter,  it  has  been  shown  that  the  Laue  technique  can  be 
reliably  applied  to  crystals  of  a  protein  as  large  as  AChE,  in  a  static  state.  The  maps  from  a  data  set 
of  24  images,  with  an  overall  completeness  of  84.5%  and  an  (unweighted)  of  14.8%,  clearly 
show  the  position  of  the  ligand  inside  the  active  site,  and  structural  changes  in  important  protein 
residues.  This  data  set  does  not  represent  the  maximum  quality  attainable  for  AChE  crystals:  use  of 
a  theta  offset  for  the  detector,  and  reduction  of  the  extent  of  X-ray  heating  of  the  crystal  could 
improve  the  useful  resolution  limit  and  the  reliability  of  the  intensities. 

To  improve  the  usefulness  of  the  Laue  method  for  AChE,  future  studies  will  concentrate  on 
synchronization  of  release  of  the  caged  compound,  and  on  coordinating  data  collection  with  either 
movement  of  choline  out  of  the  gorge  (for  the  caged  choline  compound)  or  with  carbamylation  and 
decarbamylation  (for  the  caged  CCh). 
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Structure  of  the  rhAChE/fasciculin  complex 

Knowledge  of  the  3D  structure  of  hAChE  is  of  importance  for  drug  design,  in  particular  of 
anti- Alzheimer  drugs  (102,  103),  for  development  of  improved  procedures  for  treatment  of 
intoxication  by  nerve  agents  (6),  and  for  design  of  safer  and  more  effective  insecticides. 

Vertebrate  AChE  occurs  as  an  array  of  molecular  forms  which  differ  in  the  number  of  subunits 
which  they  contain  and  in  their  quaternary  structure  (104).  Recombinant  hAChE  (rhAChE), 
expressed  in  HEK  293  cells,  is  secreted  into  the  culture  medium  as  a  mixture  of  monomers,  dimers 
and  tetramers  (105);  additional  heterogeneity  may  arise  due  to  variations  in  the  extent  of 
glycosylation  (106)  and  to  proteolytic  ‘nicking’.  Such  heterogeneity  may  severely  impede 
crystallization.  To  alleviate  this  problem,  a  mutant  of  rhAChE  was  constructed  in  which  the 
C-terminus  had  been  truncated  to  give  rise  to  a  homogeneous  monomeric  form  (107).  The  T 
subunit  expressed  was  devoid  of  the  C-terminal  hydrophobic  sequence  corresponding  to  residues 
544-583,  including  C580,  which  forms  the  interchain  disulfide  in  oligomeric  forms  (108).  These 
residues  were  replaced  by  the  pentapeptide,  ASEAP,  which  is  the  natural  C-terminal  sequence  in  the 
H  subunit  (104). 

Truncation  was  performed  by  DNA  cassette  replacement.  Specifically,  the  187  bp  DNA 
fragment  between  BssHII  and  Sail  in  pL5CA  was  replaced  by  the  synthetic  DNA  duplex.  The 
truncated  DNA  sequence  was  introduced  into  a  tripartite  expression  vector  expressing  also  the 
reporter  gene  cat  and  the  selection  marker  neo.  The  expression  vector  coding  for  the  expression  of 
monomeric  rhAChE  was  introduced  into  HEK  293  cells.  Stable  clones,  secreting  high  levels  (>  4 
U/ml/h)  of  truncated  enzyme  (ACter-rhAChE),  were  isolated  and  used  in  large-scale  production 
(109),  so  as  to  produce  amounts  of  purified  enzyme  adequate  for  crystallization.  Purification  was 
carried  out  by  affinity  chromatography  (16).  Analytical  sucrose  density  gradient  centrifugation 
showed,  as  anticipated,  that  a  homogeneous  monomeric  species  was  produced,  which  migrated 
slightly  more  slowly  than  the  WT  monomer.  SDS-PAGE  suggested  a  similar  degree  of 
glyocsylation-related  microheterogeneity  in  the  truncated  and  WT  enzymes.  The  molecular  weight 
of  the  truncated  polypeptide,  after  removal  of  sugar  side  chains  by  glycanase,  is  60kD,  as  compared 
to  64kD  for  the  WT  polypeptide.  This  difference  in  size  is  as  predicted  from  the  difference  in  chain 
length  of  the  two  forms.  Finally,  the  catalytic  characteristics  of  the  truncated  enzyme,  as  well  as  the 
inhibition  constants  for  its  inhibition  both  by  active  and  peripheral  site  ligands,  are  essentially 
identical  to  those  for  the  WT  enzyme. 

A  homologous  truncated  DNA  sequence,  in  which  site-directed  mutagenesis  afforded 
expression  of  the  E202Q  mutant  (31),  was  generated,  expressed  and  characterized  similarly  to  the 
WT  enzyme. 

Both  WT  rhAChE  and  the  E202Q  mutant  were  cocrystallized  as  a  stoichiometric  complex  with 
the  polypeptide  toxin,  fasciculin-II  (FAS-II)  (110),  from  the  venom  of  the  green  mamba,  as  had 
previously  been  done  for  both  TcAChE  (79)  and  mAChE  (78).  Crystallization  of  the  FAS-II 
complexes  of  both  WT  and  E202Q  rhAChE,  from  ammonium  sulfate  at  neutral  pH  and  19°C, 
yielded  crystals  of  irregular  prism  morphology,  which  grew  within  1-3  weeks  to  0.01-0.1  mm  in 
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size.  Recrystallization  from  water  resulted  in  larger  crystals,  up  to  0.8  mm  in  size.  The  crystals 
exhibited  rhombohedral  symmetry  of  space-group  R32,  with  unit  cell  dimensions  a=b=150A, 
c=247A,  with  a  single  copy  of  the  complex  in  the  asymmetric  unit. 

X-ray  data  were  collected  on  the  X12C  beamline  of  the  NSLS,  at  Brookhaven  National 
Laboratory,  in  both  cases  from  a  single  cryogenically  cooled  crystal.  The  overall  completeness  of 
the  30A-2.7A  data  is  99.0%  and  83.4%,  with  R,^g^gg  8.1%  and  9.1%,  for  WT  and  E202Q  AChE, 
respectively.  For  the  wild-type  enzyme,  the  current  is  22%  and  the  Rj^^^  is  29%,  and  for  the 
E202Q  structure,  the  corresponding  figures  are  21  and  24%,  respectively.  In  both  stmctures,  the 
electron  density  reveals  residues  5-490  and  497-543  of  the  rhACtiE  polypeptide  chain,  and  all  61 
residues  of  FAS-II.  The  490-497  sequence  is  homologous  to  the  external  485-490  sequence  in 
TcAChE,  which  is  not  visible  in  the  native  2ACE  stmcture,  but  can  be  seen  in  the  recently  refined 
complex  with  E2020  (1 1 1).  In  the  published  mouse  structure  (PDB  ID  IMAH),  residues  4-543  are 
all  present;  but,  judging  from  the  temperature  factors,  residues  491-499  are  poorly  determined, 
especially  at  the  3.2A  higher  resolution  limit  of  this  structure.  The  electron  density  of  the  WT 
rhAChE  structure  shows  the  location  of  -200  water  molecules  (none  were  reported  for  mAChE); 
one  A-glycosylation  site,  at  N350,  seen  also  in  the  mAChE  structure;  and  additional  density  near 
T504,  which  is  a  putative  0-glycosylation  site. 

Stmcture  solution  of  the  rhACtiE/FAS-II  complex  was  carried  out  by  molecular  replacement, 
using  the  TcAChE  /FAS-II  stmcture  (PDB  ID  IFSS,  (79)  as  a  probe  molecule  and,  subsequently, 
as  a  starting  model  for  refinement.  Diffraction  data  for  the  E202Q  mutant  were  refined  similarly, 
but  using  the  rhACtiE/FAS-II  stmcture  as  a  starting  model.  The  stmctures  were  refined  using  X- 
PLOR  (simulated  annealing  (112))  and  CNS  (simulated  annealing  coupled  with  maximum 
likelihood  technique  (113)). 

Close  inspection  of  the  rhAChE/FAS-II  stmcture  reveals  a  high  degree  of  similarity  between  it 
and  the  two  other  complexes  of  AChE  with  FAS-II,  viz-  of  TcAChE  and  mAChE  (78,  79), 
especially  in  the  active-site  region  and  within  the  aromatic  gorge.  Small  variations  can  be  found  in 
some  loops  and  in  the  insertion/deletion  regions.  Our  results  validate  use  of  TcAChE  as  a  model  for 
gaining  a  general  understanding  of  the  stmcture  of  AChE.  However,  they  also  paved  the  way  for  a 
study  of  the  subtle  differences  in  stmcture  which  may  account  for  species-dependent  differences  in 
specificity  for  substrates  and  inhibitors.  Such  subtle  differences  must  be  taken  into  account  in 
design  of  dmgs,  e.g.  for  treatment  of  Alzheimer  patients,  for  treatment  of  nerve  gas  intoxication,  and 
in  producing  insecticides  with  a  high  degree  of  selectivity  for  insect  vs.  human  AChE. 

Residue  E202,  which  is  adjacent  to  the  active-site  serine,  S203(200),  is  one  of  three  charged 
residues  at  the  bottom  of  the  active-site  gorge,  and  plays  an  important  role  in  AChE  reactivity.  Yet, 
the  3D  stmcture  of  the  E202(i99)Q  mutant  is  almost  identical  to  that  of  the  WT  enzyme.  The 
reduction  in  catalytic  efficiency  displayed  by  the  E202Q  mutant,  either  toward  charged  or 
uncharged  substrates,  and  the  marked  concomitant  decreases  in  rates  of  both  phosphylation  and 
aging  (31,  114),  was  explained  by  the  participation  of  residue  E202(i99)  in  the  chemical  process 
involving  the  catalytic  triad  residue  H447(449)  (1 15).  It  was  proposed  earlier  that  E202(i99)  plays 
a  key  role  in  maintaining  a  network  of  H-bonds  that  span  the  cross-section  of  the  active-site  gorge 
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(114).  The  role  suggested  for  this  network  is  to  stabilize  the  functional  architecture  of  the  active 
center,  and  to  ensure  proper  positioning  of  E202(/99),  thereby  stabilizing  acylation  and 
phosphylation  transition  states.  Central  to  this  array  are  the  two  water  molecules  bridging  residues 
E202(799)-E450(445)  and  E202(/99)-Y133(/3f).  Indeed,  replacement  of  any  of  these  three 
residues  reduces  catalytic  activity  to  a  similar  degree  (31).  The  two  water  molecules  appear  to  be 
conserved  in  the  E202Q  mutant.  Replacement  of  E202  by  D  leads  to  an  even  greater  decrease  in 
enzymatic  activity.  Despite  the  fact  that  the  charge  in  E202D  is  maintained,  these  data  suggest  that 
the  H-bond  network  is  dismpted  due  to  the  shorter  side  chain  of  D202.  It  will  be  of  interest  to 
investigate  whether  the  H-bond  network  involving  the  water  molecules  is  maintained  in  the  E202D 
mutant. 
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Homology  modeling  ofAChE  and  related  proteins 

The  concept  of  an  electrostatic  motif  on  the  surface  of  biological  macromolecules  as  a  definite 
topographical  pattern  of  electrostatic  potentials  in  3D  space,  provides  a  powerful  tool  for 
identification  of  functionally  important  regions  on  the  surface  of  stmcturally  related 
macromolecules  (116).  This  pattern  can  only  be  defined  in  the  context  of  molecules  of  similar  3D 
stmcture,  for  only  in  such  a  case  is  it  possible  to  unambiguously  orient  the  motif  with  respect  to 
zones  of  the  related  macromolecules  which  might  be  functionally  significant.  We  have  used  this 
approach  to  localize  a  functional  region  common  to  ChEs  and  to  a  set  of  neural  cell-adhesion 
proteins,  which  include  gliotactin  (GLI)  (117),  neurotactin  (NRT)  (118)  and  neuroligin  (NL)  (119). 
It  has  been  suggested  that  these  adhesion  proteins  are  structurally  related  to  ChEs  due  to  their  high 
sequence  similarity  (120);  however,  they  lack  the  active-site  serine. 

The  high  degree  of  identity  among  AChE  and  BChE  sequences  has  permitted  constmction  of 
homology  models  of  the  3D  stmctures  of  BChE  (33)  and  of  AChE  from  different  species  (121). 
Examination  of  the  3D  structure  of  TcAChE  showed  the  enzyme  to  be  characterized  by  a  marked 
asymmetric  spatial  distribution  of  charged  residues;  the  charges  were  shown  to  roughly  segregate 
into  a  ‘northern’  negative  hemisphere  (taking  the  mouth  of  the  gorge  as  the  north  pole)  and  a 
‘southern’  positive  one  (73).  This  electrostatic  pattern  was  shown  to  give  rise  to  a  large  dipole 
moment  which  was  calculated  to  be  -1600  Debyes,  oriented  roughly  along  the  axis  of  the  active-site 
gorge  (73, 74,  122).  Electro-optical  measures  performed  on  AChE  from  Bungarus  fasciatus  gave  a 
value  of  -1000  Debye,  thus  confirming  experimentally  the  order  of  magnitude  of  the  macrodipole 
(123). 

The  asymmetric  distribution  of  surface  potentials  was  thought  to  be  essential  for  the  fast 
catalysis  effected  by  ChEs.  This  assumption  has  been  reevaluated  in  the  light  of  site-directed 
mutagenesis  experiments  on  rhAChE  in  which  up  to  seven  charged  residues  in  the  vicinity  of  the 
entrance  to  the  active-site  gorge  were  neutralized.  The  ensuing  drastic  reduction  of  the  asymmetric 
distribution  of  surface  potentials  did  not  have  a  major  effect  on  catalysis  (124).  Brownian  dynamics 
simulations  of  bimolecular  rate  association  constants  for  ACh  and  AChE  also  failed  to  show  a 
meaningful  correlation  between  the  distribution  of  electrostatic  surface  potentials  and  catalysis 
(125).  A  further  analysis  of  the  electrostatic  properties  of  the  3D  stmctures  and  homology  models 
of  ChEs,  aimed  at  generating  a  detailed  and  quantitative  topographical  map  of  the  electrostatic 
surface  potentials,  identified  the  presence  of  an  “annular”  electrostatic  motif  of  negative  surface 
potential  around  the  entrance  of  the  active-site  gorge  (121).  The  role  of  surface  potentials  in 
catalysis  is  not  clear-cut.  It  has  been  shown  to  play  only  a  minor  role  in  the  case  of  isozymes  of 
trypsin  differing  by  as  much  as  12.5  charge  units.  The  catalytic  efficiency  of  these  isozymes  was 
shown  to  be  influenced  not  by  electrostatic  surface  potentials,  but  by  a  few  charged  amino  acids 
localized  close  to  the  active  site  (126).  These  considerations  drove  us  to  investigate  whether  the  role 
of  surface  potentials  in  ChE  might  be  related  to  roles  other  than  catalysis. 

ChEs  possess  a  3D  fold  different  from  those  of  “classical”  serine  hydrolases,  such  as 
chymotrypsin  and  subtilisin.  This  fold  is  called  the  a/p  hydrolase  fold  (127),  and  has  been  shown 
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to  characterize  a  family  of  hydrolytic  enzymes  of  widely  different  phylogenetic  origin  and  substrate 
specificity.  Some  of  these  enzymes,  like  the  fungal  triacylglycerol  lipases  from  Geotrichum 
candidum  and  Candida  rugosa  (128),  show  substantial  sequence  similarity  with  ChEs.  Other 
members  of  the  family,  like  haloalkane  halogenase  from  Xenobacter  autotrophicus  and 
dienelactone  hydrolase  from  Pseudomonas  sp.,  do  not  show  significant  sequence  similarity  to  the 
ChEs;  but  inspection  of  their  3D  stmctures  shows  that  their  polypeptide  chain  has  the  stmctural 
core  typical  of  the  a/p  hydrolase  fold  (127).  A  systematic  search  of  sequence  databases  has 
retrieved  a  large  number  of  proteins  with  a  high  degree  of  sequence  similarity  to  a/p  hydrolases 
(120).  Inspection  of  the  sequences  reveals  that  amongst  an  overwhelming  majority  of  enzymes, 
there  are  also  a  few  proteins  involved  in  cell  adhesion;  these  contain  a  ChE-like  domain  whose 
function  is  most  probably  independent  of  hydrolase  activity  since  the  active-site  serine  is  mutated  to 
a  glycine.  Amongst  them  are  the  three  proteins  which  we  chose  to  investigate,  because  of  their 
presumed  involvement  in  cell  adhesion  and  signaling  in  the  nervous  system. 

GLI  is  a  transmembrane  protein  transiently  expressed  on  peripheral  glia  of  Drosophila 
melanogaster  ercAixyos  and  required  for  formation  of  the  peripheral  blood-nerve  barrier  (117).  The 
C-terminal  extracellular  domain  of  this  protein  shows  substantial  sequence  similarity  with  ChEs. 
Such  similarity  can  also  be  found  in  the  extracellular  domain  of  NRT,  a  transmembrane 
glycoprotein  which  is  expressed  in  neuronal  and  epithelial  tissues  during  embryonic  and  larval 
stages  (118,  129).  Sequence  similarities  suggest  that  neuroligins  (NLs),  a  family  of  mammalian 
neuronal  cell  surface  protein  that  act  as  ligands  for  P-neurexins  (119),  also  contain  an  extracellular 
ChE-like  domain.  NLs  and  neurexins  mediate  interactions  between  neurons  and  contribute  to  the 
specific  organization  of  synapses.  Both  the  Drosophila  adhesion  proteins  and  NL  are  involved  in 
heterophilic  cell-cell  interactions,  suggesting  that  their  ChE-like  domains  possess  a  recognition 
function  that  has  been  conserved  from  invertebrates  to  vertebrates.  It  is  worth  noticing  that  both  GLI 
and  NRT  display  susbstantially  higher  sequence  similarity  with  TcAChE  than  with  the  DmAChE 
enzyme  (Table  3),  and  that  the  targets  of  their  adhesive  functions  have  not  yet  been  identified. 

The  sequence  similarity  between  ChEs,  NL,  GLI,  and  NRT  prompted  us  to  build  homology 
models  of  their  extracellular  domains,  both  to  confirm  the  presence  of  shared  structural  motifs  and 
to  investigate  the  possibility  that  they  might  display  similar  electrostatic  characteristics.  Our 
underlying  assumption  is  that  the  ChE-like  domains  of  these  adhesion  proteins  are  independently 
folded,  and  that  neither  their  conformation  nor  their  electrostatic  properties  will  be  greatly  affected 
by  the  association  with  the  non-ChE  Uke  regions  of  GLI,  NL  and  NRT,  which  are  mostly 
transmembrane  and  intracellular  (Fig.  10). 

The  homology  models  were  built  with  the  automated  knowledge-based  building  tool  Swiss- 
Model  (130),  and  validated  as  described  previously  (121).  The  same  modelling  procedure  was 
applied  as  a  control  to  the  neutral  lipase  from  G.  candidum  (GcLip)  (131),  whose  3D  structure  is 
known,  to  check  for  the  bias  in  modelling  introduced  by  the  templates.  The  modelling  of  GcLip 
gave  a  structure  very  close  to  the  actual  3D  stmcture  solved  by  X-ray  crystallography  (a-carbon 
RMS  deviation  ~1.1A  from  the  3D  X-ray  stmcture  of  GcLip).  Since  all  the  proteins  involved  in  this 
study  possess  about  the  same  degree  of  sequence  identity  as  GcLip  does  with  TcAChE  (see  Table 
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3),  we  can  reasonably  assume  our  models  to  be  close  approximations  of  the  actual  structures  of 
these  domains. 

The  electrostatic  surface  potentials  of  TcAChE  (PDB  entry  2ACE),  GLI,  NL,  and  NRT  were 
calculated  solving  the  Poisson-Boltzmann  equation  by  the  finite  difference  method  (FDPB)  (132), 
using  the  DelPhi  algorithm  (133)  as  implemented  in  the  GRASP  program  (134).  As  a  control,  we 
examined  the  distribution  of  surface  potentials  for  other  members  of  the  a/p  hydrolase  family  to 
check  whether  the  “annular”  electrostatic  motif  was  a  characteristic  shared  by  most  members  of 
this  family.  We  chose  three  enzymes  with  a  sequence  identity  to  ChEs  similar  to  that  of  NL,  NRT 
and  GLI.  The  three  enzymes  selected  were  cholesterol  esterase  from  Candida  rugosa  (135)  (CrCE, 
PDB  entry  ICLE),  GcLip  (PDB  entry  ITHG)  and  hile-salt-stimulated  lipase  from  rat  (rBSSL) 
(136).  In  this  last  case  the  3D  stmcture  is  not  yet  available,  and  a  homology  model  was  built 
following  the  same  procedure  used  for  GLI,  NRT  and  NL.  We  also  calculated  the  surface  potentials 
for  hBChE,  for  hAChE  and  for  the  mutant  of  rhAChE  in  which  seven  negatively  charged  residues 
near  the  entrance  of  the  gorge  had  been  neutralized  (7-hAChE).  The  structures  of  7-hAChE,  NL, 
NRT  GLI,  GcLip,  CrCE  and  rBSSL  were  superimposed  on  that  of  TcAChE  by  a  least-squares 
fitting  procedure  implemented  in  the  program  LSQMAN  (137),  and  placed  in  a  common  orientation 
in  which  the  axis  of  the  active-site  gorge  of  TcAChE  (as  defined  by  Antosiewicz  et  al.  (122))  is 
aligned  with  the  cartesian  Z-axis. 

Our  most  significant  result  is  shown  in  Fig.  11.  It  is  immediately  apparent  that  the  electrostatic 
surface  potential  contours  of  TcAChE,  GLI,  NRT  and  NL  are  strikingly  similar.  All  three  proteins 
show  a  characteristic  “annular”  electrostatic  motif  of  negative  potential  around  the  zone 
homologous  to  the  active-site  gorge  of  the  ChEs.  In  contrast,  GcLip,  rBSSL,  and  CrCE  appear  to 
display  a  different  topographical  distribution  of  surface  potentials  (Fig.  1  la).  In  the  case  of  the  7- 
hAChE  mutant,  neutralization  of  seven  negatively  charged  residues  in  the  “annular”  region 
virtually  erases  the  electrostatic  motif,  while  not  affecting  the  negative  potential  gradient  within  the 
active-site  gorge.  This  distribution  of  electrostatic  surface  potentials  is  very  similar  to  that  of  hBChE 
(Fig.  11b).  This  finding  correlates  weU  with  the  contention  that  the  potential  gradient  inside  the 
active-site  gorge  makes  a  more  important  contribution  to  ChE  catalysis  than  the  surface  potentials 
around  the  rim  of  active-site  gorge  (124,  138). 

In  order  to  provide  a  quantitative  measure  of  the  electrostatic  motif  of  the  molecules  studied,  we 
calculated  the  average  surface  potential  over  a  20A  wide  region  of  accessible  surface  around  the 
zone  homologous  to  the  entrance  of  the  active-site  gorge  of  TcAChE  (which  constitutes  the 

o 

common  Z  axis).  This  region  was  divided  into  lA  wide  concentric  sections,  and  the  average 
potentials  in  the  sections  were  plotted  and  correlated  by  a  linear  regression  fit  with  the  results 
obtained  for  TcAChE  (see  Table  4,  Fig.  12).  GLI,  NRT  and  NL  are  all  characterized  by  a  high 
correlation  coefficient  and  an  average  potential  in  the  “annular”  region  very  similar  to  that  of 
TcAChE;  in  addition,  their  potential  gradient  closely  follows  that  of  TcAChE.  Of  the  three  lipases, 
only  CrCE  shows  a  high  correlation  coefficient  and  an  average  potential  similar  to  AChE,  but  its 
RMS  deviation  from  the  ChE-like  fold  is  the  largest  amongst  all  the  proteins  examined,  and  its 
potential  gradient  is  substantially  different  from  that  of  TcAChE.  It  is  of  interest  that  BChE  displays 
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a  lower  average  potential  in  the  “annular”  region  and  a  very  poor  correlation  coefficient,  in  addition 
to  a  potential  gradient  similar  to  that  of  the  lipases. 

The  biological  significance  of  the  electrostatic  surface  potential  of  ChEs  has  been  the  subject  of 
speculations,  and  it  has  been  proposed  that  it  might  have  a  non-catalytic  function  (124).There  is 
mounting  evidence  that  ChEs  may  be  involved  in  functions  distinct  from  their  catalytic  activity. 
During  development  of  the  nervous  system  of  vertebrates,  expression  of  ChEs  is  associated  with 
periods  of  enhanced  neurite  outgrowth  (139,  140),  suggesting  that  they  might  serve  as  neuronal 
growth  factors  or  cell-adhesion  molecules.  It  has  been  speculated  that  the  site  responsible  for  the 
non-catalytic  interaction  of  ChEs  is  localized  close  to  the  entrance  of  the  active-site  gorge,  possibly 
overlapping  with  what  has  been  termed  the  “peripheral”  anionic  site  (PAS)  of  AChE,  a  possible 
locus  of  regulation  of  enzymatic  activity  (27,  30,  141),  which  is  located  in  the  region  of  the 
“annular”  electrostatic  motif.  In  a  recent  study,  AChE  was  shown  to  promote  neurite  regeneration 
in  cultured  adult  neurons  of  the  mollusk  Ap/y^/a  (142).  The  growth-promoting  action  of  AChE  was 
shown  to  be  independent  of  its  catalytic  function.  Inhibition  of  the  active  site  did  not  influence  the 
neurotrophic  action  of  AChE,  but  occupancy  of  its  PAS  suppressed  neurite  regeneration  (142). 
There  is  also  evidence  that  the  PAS  of  AChE  can  enhance  the  rate  of  formation  of  amyloid  fibrils,  a 
major  component  of  the  senile  plaques  found  in  the  brains  of  Alzheimer  patients  (143).  Strong 
evidence  for  a  common  recognition  mechanism  in  cell-adhesion  between  ChEs,  GLI  and  NLs 
comes  from  a  recent  study,  in  which  chimeric  constmcts  of  the  C-terminal  transmembrane  domain 
of  NRT,  fused  to  the  extracellular  ChE-like  domains  of  DmAChE  or  TcAChE,  were  shown  to  be 
endowed  with  the  same  heterophilic  adhesion  properties  as  wild  type  NRT  (144). 

Electrostatic  complementarity  between  interacting  proteins  has  been  found  to  be  one  of  the 
major  driving  forces  for  complex  formation  (145).  Mutagenesis  and  electrostatic  screening  studies 
of  the  association  rates  between  bamase,  a  ribonuclease,  and  barstar,  its  inhibitor,  have  shown  the 
limiting  step  to  be  the  formation  of  a  binding  transition  state,  in  which  the  two  proteins  form  a  low- 
affinity  non-specific  complex  held  together  by  long-range  electrostatic  interactions  (146).  We 
speculate  that  a  possible  mechanism  effected  by  the  “annular”  motif  might  be  the  facilitation  of 
target  recognition  via  low  specificity  associations  mediated  by  long-range  electrostatic  interactions 
with  a  pool  of  ligands  displaying  a  complementary  motif,  followed  by  a  subsequent  docking  step 
with  the  unique  functional  ligand(s),  to  yield  the  final  high  affinity  complex.  Some  support  for  our 
hypothesis  comes  from  the  mode  of  binding  of  AChE  to  one  of  its  strongest  inhibitors,  fasciculin  II 
(FAS),  a  three-fingered  polypeptide  toxin  from  the  venom  of  the  green  mamba  (110).  These 
polypeptide  snake  venom  toxins  constitute  a  family  of  6-8  kDa  proteins,  of  which  a-bungarotoxin, 
a  potent  antagonist  for  the  nicotinic  ACh  receptor,  is  the  best  known  member. 

FAS  has  been  shown  to  bind  in  a  “cork-and-bottle  fashion”  to  the  entrance  of  the  active-site 
gorge  of  AChE,  which  is  situated  in  the  negatively  charged  “northern”  hemisphere  (78,  79).  The 
surface  complementarity  between  enzyme  and  inhibitor  is  very  large;  about  2000A2  of  accessible 
surface  are  buried  upon  binding,  including  the  region  of  the  “annular”  motif  In  addition  to  this 
large  surface  complementarity,  the  AChE-FAS  complex  displays  significant  electrostatic 
complementarity.  Analysis  of  the  topography  of  surface  potentials  of  FAS  shows  that  it  also 
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possess  a  “bipolar”  motif,  with  a  preponderance  of  positive  surface  which  is  absent  in  other  three- 
fingered  snake  toxins  (110).  This  positive  surface  characterizes  the  zone  of  interaction  between 
FAS  and  AChE  (79). 

Analysis  of  the  sequence  alignment  of  the  proteins  which  we  investigated  shows  no  sequence 
motif  shared  exclusively  by  AChE,  GLI,  NRT  and  NL.  Thus  the  conjunction  of  the  topographical 
distribution  of  surface  potentials  with  the  topological  arrangement  of  the  polypeptide  chain  seems 
the  most  important  common  stmctural  factor  conserved  between  these  proteins.  It  is  thus  plausible 
that  the  shared  electrostatic  motif  indicates  a  common  recognition  mechanism  or  even  a  common 
ligand.  We  hypothesize  that  a  chimeric  construct  of  the  cytoplasmic  domain  of  NRT  with  the 
mutant  form  of  AChE,  in  which  seven  negative  residues  have  been  neutralized  to  abolish  the 
electrostatic  motif,  will  not  display  the  same  adhesive  properties  as  a  construct  built  with  wild-type 
AChE.  We  further  predict  that  constmcts  with  BChE  and  CrCE  will  not  be  able  to  promote 
heterogeneous  cell  adhesion,  the  former  because  it  lacks  the  “annular”  motif,  and  the  latter  because 
of  a  large  deviation  from  the  common  ChE-like  fold. 

The  above  considerations  lead  us  to  conclude  that  GLI,  NRT,  NLs  and  AChE  are  the  first 
members  of  a  class  of  adhesion  proteins  that,  because  of  their  common  electrostatic  and  stmctural 
motif,  we  have  decided  to  name  “electrotactins”. 
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Studies  on  human  paraoxonase 

Human  paraoxonase/arylesterase  (PONl)  is  an  esterase  with  broad  substrate  specificity.  It 
catalyzes  the  hydrolysis  of  paraoxon,  and  a  number  of  other  OPs,  and  of  esters  such  as 
phenylacetate.  OPs  are  widely  used  as  insecticides  and  nerve  gases,  and  are  lethal  for  many  species, 
including  humans,  due  to  their  inhibition  of  the  synaptic  enzyme,  AChE,  as  well  as  of  other 
esterases  and  proteases.  It  has  been  shown  that  serum  paraoxonase  plays  a  significant  role  in  the 
detoxification  of  a  variety  of  OPs  in  mammals.  Furthermore,  correlations  have  been  found  between 
the  levels  of  semm  paraoxonase  and  systemic  diseases,  such  as  systemic  amyloidosis  (see 
BACKGROUND).  The  enzyme  does  not  correspond  to  any  of  the  protein  sequences  in  the  Swiss- 
Prot  data  base  or  in  the  protein  sequences  derived  from  the  GENEBANK  database;  thus  it  is 
difficult  to  infer  any  structural  information  from  proteins  found  in  the  Brookhaven  Protein  Data 
Bank  (PDB)  (147).  Although  PONl  may  be  involved  in  breakdown  of  oxidized  lipids 
(BACKGROUND),  its  physiological  role  and  natural  substrate  remain  to  be  clearly  defined. 
Structural  data  may  shed  light  on  this  issue,  and  may  also  offer  a  new  structural  motif 

Expression  of  PONl  has  been  achieved,  both  in  human  293  kidney  cells  and  in  CHO-Kl  cells. 
These  expression  systems  were  utilized,  in  conjunction  with  site-directed  mutagenesis,  to  obtain 
evidence  for  a  role  in  the  catalytic  action  of  PON  1  of  various  candidate  amino  acids,  and  in  an 
attempt  to  increase  its  catalytic  efficacy  (148,  149).  In  order  to  obtain  the  large  quantities  of  pure 
protein  necessary  for  X-ray  crystallographic  studies,  both  of  wild  type  PONl  and  of  pertinent 
mutant  forms  of  the  enzyme,  we  have  chosen  the  Baculovims  (Bac)  expression  system  in  insect 
cells.  This  expression  system  has  proved  itself  in  recent  years  to  be  an  extremely  versatile  and  cost 
effective  tool  for  large  scale  production  of  many  eukaryotic  proteins,  with  yields  of  as  much  as 
50mg/l  of  culture  medium  (150).  A  scan  of  the  PDB  retrieved  almost  100  structures  of  protein 
which  had  been  expressed  using  the  Bac  expression  system  fhttp://www.pdb.bnl.gov').  The  Bac 
system  has  also  been  successfully  used  in  expression  of  both  wild-type  and  selenomethionyl 
variants  of  human  choriogonadotropin  (151).  The  production  of  selenomethionyl  variants  of  a 
particular  protein  is  an  extremely  powerful  tool  which  permits  determination  of  a  structure  by 
multiple  anomalous  diffraction  (MAD)  techniques  (152).  This  technique  requires  only  a  single 
crystal,  with  the  phase  problem  being  solved  by  varying  the  wavelength  during  data  collection.  This 
is  in  contrast  to  multiple  isomorphous  replacement  methods,  which  require  crystals  that  are 
isomorphous  to  the  native  ones  with  heavy  atoms  soaked  in. 

The  Bac-based  system  is  a  eukaryotic  expression  system;  it  thus  uses  many  of  the  protein 
modification,  processing,  and  transport  systems  present  in  higher  eukaryotic  cells.  Bac  is  a  large, 
enveloped  double-stranded  DNA  vims  that  infects  arthropods.  The  viral  genome  is  large  (130  kbp); 
hence  it  can  accommodate  large  segments  of  foreign  DNA.  Vimses  enter  the  cell  by  adsorptive 
endocytosis  and  move  to  the  nucleus,  where  their  DNA  is  released.  Two  types  of  viral  progeny  are 
produced:  extracellular  vims  is  released  from  the  cell  by  budding,  starting  ~12  hr  postinfection,  and 
is  produced  at  a  logarithmic  rate  until  20  hr  postinfection,  after  which  production  falls  off 
Polyhedra-derived  vims,  appears  in  the  nucleus  at  ~18  br  postinfection,  and  continues  to  accumulate 
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as  late  as  72  hr  postinfection,  or  until  the  cells  lyse.  Occluded  viral  particles  are  embedded  in 
proteinaceous  viral  occlusions  called  polyhedra  within  the  nuclei  of  infected  cells.  The  polyhedrin 
protein  (29  kDa)  is  the  major  protein  component  of  the  occlusion  bodies. 

The  Bac  expression  system  takes  advantage  of  several  properties  of  the  polyhedrin  protein:  (1) 
It  is  expressed  at  very  high  levels  in  infected  cells,  constituting  more  than  half  of  the  total  cellular 
protein  late  in  the  infectious  cycle  ;  (2)  It  is  not  essential  for  infection  or  replication  of  the  vims;  (3) 
Vimses  lacking  the  polyhedrin  gene  have  a  plaque  morphology  distinct  from  that  of  vimses 
containing  the  gene.  Thus  recombinant  Bac  is  engineered  to  express  the  desired  protein  by 
replacing  the  polyhedrin  gene  with  a  foreign  gene  coding  for  this  protein  through  homologous 
recombination.  This  is  achieved  by  contransfecting  a  linearized  Bac  genome  with  a  transfer  vector 
containing  the  gene  of  interest.  This  transfer  vector,  besides  the  elements  needed  for  amplification  in 
E.  coli,  contains  a  multiple  cloning  site  for  the  insertion  of  the  target  gene,  a  promoter  for  the 
polyhedrin  gene  and  a  large  tract  of  Bac  genome  flanking  the  cloning  region  to  facilitate 
homologous  recombination.  Subsequently,  recombinant  vimses  in  which  expression  of  polyhedrin 
has  been  replaced  by  the  target  protein  can  be  selected  by  several  screening  methods.  Many 
commercial  improvements  to  this  basic  protocol,  that  report  increased  efficiency  in  the  produttion  of 
recombinant  vims,  are  available  on  the  market,  thus  minimizing  the  time  needed  for  production  of  a 
suitable  vims. 

We  have  received,  from  the  laboratory  of  Bert  La  Du  (Department  of  Pharmacology,  University 
of  Michigan),  with  whom  we  are  collaborating,  a  variant  of  the  PONl  gene  in  which  residues  H19 
and  Q20  have  both  been  mutated  to  Ala.  This  constmct  permits  secretion  of  PONl  into  the 
medium,  by  enabling  the  cell  machinery  to  cleave  the  N-terminal  signal  peptide  (153).  While  the 
resulting  product  displays  reduced  activity,  its  solubility  properties  are  greatly  enhanced  and  we  feel 
that  it  will  constitute  an  ideal  candidate  for  stmctural  studies.  Production  of  recombinant  PON  1  will 
be  monitored  using  antibodies  supplied  by  the  La  Du  laboratory. 

In  preparatory  experiments,  performed  in  our  laboratory,  the  PON  1  clone  received  from  the  La 
Du  laboratory  has  been  amplified  via  PCR,  in  a  BlueScript®  shuttle  vector,  with  primers  designed  to 
contain  the  appropriate  restriction  sites  for  unidirectional  cloning  in  the  transfer  vector.  A  region 
coding  for  six  histidine  residues  (His-tag)  was  included  in  the  3'  primer,  in  order  to  permit 
purification  of  the  expressed  protein  on  a  nickel  column  (154).  The  product  was  sequenced  so  as  to 
confirm  the  integrity  of  the  coding  sequence,  and  to  check  whether  PCR  had  introduced  any  errors. 
The  gene  fragment  was  then  isolated  by  digestion  with  the  appropriate  restriction  enzymes,  purified 
and  cloned  in  an  appropriate  transfer  vector  for  the  Bac  expression  system,  viz.  AcNPV  (155).  The 
initial  attempts  to  express  PON  1  in  Bac  will  be  performed  in  the  next  few  weeks. 
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TABLES 


Table  1.  Monochromatic  diffraction  data  statistics  for  TcAChE-edrophonium  complex 


No.  of  frames 

66 

Oscillation  angle 

0.75° 

Cell  parameters  (P3,21) 

a=b=112.41A,  c=136.55A 

Mosaicity 

0.29° 

No.  of  observations 

275,832 

Unique  reflections 

38,611 

Overall  redundancy  (highest  resolution  shell*) 

3.1  (3.1) 

(highest  resolution  shell) 

8.9%  (42.0%) 

Average  I  /  average  a  (highest  resolution  shell) 

14.0  (2.1) 

Completeness  (highest  resolution  shell*) 

97.7%  (99.6%) 

*highest  resolution  shell:  2.46  -  2.40A 
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Table  2.  Laue  X-ray  diffraction  data  statistics  for  TcAChE-edrophonium  complex 


No.  of  frames 

24 

Exposure  time 

1  ms 

Cell  parameters  (P3,21) 

a=b=113.00A,  c=136.80A 

Singles 

No.  of  observations* 

83,171 

Unique  reflections 

19,613 

Overall  redundancy 

4.2 

■^sym 

Unweighted 

14.8% 

Weighted 

10.8% 

Singles  and  multiples  combined 

No.  of  deconvoluted  multiples 

5,464 

Total  No.  of  unique  reflections 

21,122 

Resolution  (A) 

Completeness!  (%) 

Rsy™  (%) 

Shell  Cumulative 

oo-  11.2 

42.0  42.0 

14.7 

11.2-8.4 

77.6  62.7 

13.0 

8.4  -  5.6 

91.2  82.8 

12.8 

5.6  -  4.0 

95.1  90.6 

13.5 

4.0  -  3.4 

95.3  92.4 

15.2 

3.4  -  3.0 

87.3  90.4 

21.3 

3.0  -  2.8 

61.3  84.5 

25.0 

*  After  rejection  of  weak  spots  and  those  with  negative  intensity  (I/a(I)  <  0.1),  high  background,  bad 
profile  fitting,  and/or  bad  merging 

t  Completeness  by  resolution  shell  for  singles  and  multiples  combined. 
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Table  3.  Sequence  identity  and  similarity  between  electrotactins,  lipases  &  ChEs 


Protein 

Percent  identity 

Percent  similarity 

to  TcAChE 

to  TcAChE 

(DmAChE) 

(DmAChE) 

Cholinesterases 

TcAChE 

100.0 

100.0 

tiACtiE 

57.6 

85.0 

mAChE 

59.0 

85.3 

DmAChE 

36.3 

65.8 

7-hAChE 

57.6 

85.0 

hBChE 

52.4 

83.1 

Adhesion  proteins 

NRT  (Drosophila) 

31.1  (25.1) 

57.7  (50.4) 

GLI  (Drosophila) 

30.6  (23.5) 

61.0(55.7) 

NL  (mouse) 

32.3 

60.1 

Lipases 

GcLip 

29.4 

56.8 

rBSSL 

31.2 

61.4 

CrCE 

31.9 

56.8 

Percent  sequence  identity  and  similarity  between  electrotactins,  lipases  and  ChEs,  as  given  by  the 
BestFit  algorithm  from  the  GCG  package  (version  9,  Genetics  Computer  Group,  575  Science  Drive, 
Madison,  WI 53711,  US  A)  .The  values  are  calculated  with  respect  to  the  sequence  of  rcACtiE. 
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Table  4.  Structural  deviation  between  electrotactins,  lipases  and  ChEs 


Protein 

RMS  deviation  (A)  of 
Ca  atoms  relative  to 
TcAChE  (%  of 

matched  residues) 

Average  surface 
potentials  (kT/e) 
in  the  “annular” 
motif  region 

Correlation  coefficients 
for  surface  potentials  in 
concentric  sections  of 

the  “annular”  motif 
region 

Cholinesterases 

TcAChE 

0.00  (100) 

-2.272 

1.00 

hAChE 

0.48  (97.9) 

-1.139 

0.89 

rnACtiE 

0.31  (97.1) 

-1.409 

0.80 

DmAChE 

0.77  (98.1) 

-2.358 

0.86 

7-hAChE 

0.48  (97.9) 

-0.052 

0.90 

hBChE 

0.42  (94.9) 

-0.724 

0.66 

Adhesion  proteins 

Nrt  (Drosophila) 

1.19(94.7) 

-2.132 

0.83 

Gli  (Drosophila) 

1.12(86.5) 

-1.194 

0.94 

NL  (mouse) 

1.10  (84.6) 

-1.332 

0.89 

Lipases 

GcLip 

1.46  (64.4) 

-0.551 

-0.45 

rBSSL 

0.85  (93.0) 

-0.020 

0.59 

CrCE 

1.54  (71.5) 

-1.571 

0.89 

RMS  structural  deviation  (A)  of  the  Ca  atoms,  average  surface  potentials  (kT/e)  and  correlation 
coefficients  for  electrostatic  surface  potentials  between  electrotactins,  lipases  and  ChEs.  The  values 
are  calculated  with  respect  to  the  3D  structure  and  electrostatic  properties  of  TcAChE. 
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FTGURRS 


Fig  1 .  Difference  Fourier  map  of  the  'aged'  soman-TcAChE  conjugate 


- 


Fig.2.  Active  site  of  MeP-TcAChE,  obtained  by  reaction  with  sarin 

Note  the  proximity  of  four  H-bond  donors  (dashed  lines)  to  the  anionic  phosphonyl  oxygen  atoms:  viz  the  backbone 
amide  nitrogens  of  A201,  G118,  andG119,  andN£2  of  H440.  The  OP  methyl  carbon  is  within  non-bonded  contact 
distances  (dotted lines)  of  F288  (3,6A  to  Cel)  and  F290  (3.6A  to  Cel)  in  the  acyl  pocket.  For  reference,  the  Ca- 
traces  of  native-  and  OP-enzyme  are  overlaid.  The  active  site  shown  is  essentially  identical  to  that  obtained  with 
soman  (not  shown),  except  for  minor  differences  in  inter-atomic  distances. 


Fig.3.  Active  site  of  MiPrP-TcAChE  obtained  by  reaction  with  DFP 

As  with  aged  phosphonylated  AChE,  the  stability  of  M/PrP-AChE  apparently  stems  from  four  potential  H-bond 
donors  (dashed  lines)  to  the  anionic  phosphoryl  oxygen  atoms.  However,  unlike  for  soman  or  sarin,  the  reaction  with 
DFP  results  in  a  distortion  of  the  main  chain  in  the  acyl  pocket  of  AChE.  F288  and  F290  move  significantly  from 
their  positions  in  the  native  enzyme  (shown  with  thin  dark  lines  for  reference).  The  Ca-  traces  of  native-  and  MiPrP- 
TcAChE  are  overlaid. 
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Edrophonium 


Fig.  4.  Structure  of  edrophonium 
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Fig.  5.  Distribution  of  single  reflections:  resolution  versus  wavelength 

The  Bragg-angle  cut-off  is  clearly  seen  at  higher  wavelengths  and  higher  resolution. 
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Fig.  6.  Normalization  curve  for  ID9,  ESRF,  deduced  from  the  Laue  data  set  for 
TcAChE 


(mono) 


Fig.  7.  Correlation  between  Laue  and  monochromatic  structure-factor  amplitudes 

The  correlation  is  91.6%,  based  on  19,725  unique  reflections  that  were  observed  in  both  data  sets. 
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Fig.  8.  Difference  electron  density  maps  for  the  EDR  in  the  EDR/TcAChE 
complex 


8a)  All  monochromatic  data  to  2.4A.  The  ligand  is  emphasized  with  larger  balls  and  sticks.  maps  at  3.5 

CT  cut-off  are  shown  in  green;  only  electron  density  around  the  ligand  is  shown. 


8b)  Selected  monochromatic  reflections  which  are  also  present  in  the  Laue  data  set  (2.8A).  Representation  as  in  Fig. 
8a. 
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Fig.  9,  Conformational  change  in  S200  produced  upon  binding  of  edrophonium 


9a)  Difference  electron  density  maps  from  the  Lane  data  directly  after  rigid-body  refinement 

pbs  pale  drawn  in  the  vicinity  of  Ser200,  are  contoured  at  2.7  and  -2.7  o  cut-off,  and  shown  in  green  and  red, 
respectively.  Protein  residues  in  the  native  structure  are  shown  in  black. 


9b)  Stick  model  of  (i)  the  native  starting  model  (shown  in  black);  (ii)  refined  model  using  monochromatic  data 
(shown  in  green);  (iii)  refined  model  using  Laue  data  (shown  in  gold). 
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Fig.lO.  Schematic  representations  of  domain  organizations  of  GLI,  NRT  &  NL 
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Fig  J 1  .GRASP  image  of  electrostatic  surface  potentials  in  elect  rotactins  & 
lipases 


lla)  (a)  TcAChE;  (b)  NRT;  (c)  NL;  (d)  GLI;  (e)  CrCE;  (0  rBSSL;  (g)  GcLip.  The  presence  of  the  “annular"’ 
electrostatic  motif  in  the  adhesion  molecules  and  in  AChE  should  be  noted.  The  values  of  surface  potentials  are 
expressed  as  a  spectrum  ranging  from  -\-0.25kTle  (deep  blue)  through  0  kT/e  (white)  to  -0.25  kT/e  (deep  red).  At  a 
temperature  of  298. 15  K,  kT/e  =  25.7  mV.  The  surface  potentials  are  calculated  for  a  salt  concentration  of  0. 145  M, 
with  a  cubic  grid  of  linear  dimensions  of  65A,  a  protein  dielectric  constant  of  4  and  a  solvent  dielectric  constant  of 
78. 
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lib)  (a)  hAChE;  (b)  TcAChE;  (c)  7-hAChE;  (d)  hBChE,  showing  the  absence  of  the  “annular”  motif  in  hBChE  and 
its  neutralization  in  the  7-hAChE  mutant. 
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Fig.l2. Plots  of  average  surface  potentials  measured  in  concentric  “annular’* 
sections. 

The  plots  show  the  change  in  average  potential  perpendicular  to  the  gorge  axis,  in  the  “annular”  motif  region,  for 
hBChE,  CrCE,  GcUp,  BSSL,  NL,  NRT,  GLI  and  TcAChE.  The  origin  is  placed  on  the  axis  of  the  active-site  gorge 
of  TcAChE  (122),  and  the  potential  values  are  expressed  in  units  of  kite,  where  k  is  Boltzmann’s  constant,  T  the 
temperature  in  K,  and  e  the  electronic  charge.  At  a  temperature  of  298.15  K,  kTle  =  25.7  mV. 
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